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ABSTRACT: Independent Modal Space Control ( JMSC ) is a central 

scheme that decouples the space structure into n independent 
second-order subsystems according to n controlled modes and 
controls each mode independently. It is well-knc*n that the BtSC 
eliminates control and observation spillover caused when the 
conventional coupled modal control scheme is employed . The 
independent control of each mode requires that the number of 
actuators be equal to the number of modelled modes, which is 
very high for a faithful modelling of large space structures. In 
this paper, me propose a control scheme that allows one to use a 
reduced number of actuators to control all modelled modes 
subcptimally. In particular, the method of generalized inverse 
matrices is employed to ijqilement the actuators such that the 
eigenvalues of the closed-loop system are as closed as possible 
to those specified by the optimal DISC. Computer simulation of 
the proposed control scheme on a simply supported beam is given. 

l. iNnmjcncH 

The development of the space shuttle has opened the possibility 
of constructing very large structures in space for space 
explorations. Two control problems for LSS are attitude control 
and shape control. Complex missions impose many stringent 
requirements on shape and attitude of the LSS, which lead the 
control researchers to the concept of distributed active control 
where several actuators and sensors are placed on the structure 
to in order to optimize its performance and behavior. There has 
been a considerable interest in the area of active control of 
large space structures (LSS) [1J~[13]. A number of control 
schemes were studied, but they represent one form or another of 
modal control [6]. Two main modal control schemes are the 
coupled modal control and the Independent Nodal Space Control 
(DOC) . The former uses an active controller that consists of a 
state estimator and a state feedback; the latter decouples the 
LSS into n independent subsystems according to n controlled modes 
and controls each mode independently by means of a modal filter 
[ 5 ] and an optimal controller. Coupled modal control causes 
control and observation spillover that together can destabilize 
the LSS [10]. DISC does not have the spillover problem since 
each mode is controlled independently. However in order to 
implement the DISC the number of actuators is required to be 
equal to the mmber of controlled modes which is usually very 
huge for a faithful modelling of the LSS. This fact presents a 
fundamental limitation of DISC since the required number of 
actuators is unrealizable. The main objective of this paper is to 
implement the DISC with a milder requirement of the actuator 
mmber. In other words, we will develop a control scheme that 
uses a reduced mmber of actuators to control all modelled modes 
in such a way that the modes of the closed-loop system are as 
closed as possible to the optimal nodes specified by the DCSC 
scheme. In particular, the method of generalized inverse 
matrices is employed for the implementation of DISC. 

Matrix notations used in this paper is given below: 


(1) Block diagCM}, M 2 , ...,1^) 


(2) 0 ro * « n null matrix 

(3) 1^ *um identity matrix 


Ml 0 
0 Mj 


0 0 


2 . sutojer w hswbr Km swx ocnidcl 


The description of a large flexible space structure is given ky 
the followim partial dif f erental differential equations [3] : 

m(p) aW.o/a t 2 + ua(p.t) - f(p,t) (i) 

that must be satisfied at every point P of the domain D, where 
u(P,t) is ti* displacement of Point P, La linear differential 
self-adjoint operator of order 2p, expressing the system 
stiffness, MIP) the distributed mass, and F(P,t) the distributed 
control force. The displacement u(P,t) is subject to the 
boundary conditions: 

TiuCP't) - 0; i » 1,2, ...,p (2) 

where Tj, i*l,2,...,p are linear differential operator of order 
ranging Iron 0 to (2p-l). 

The associated eigenvalue problem is formulated by: 

L* r (P) - XM(P) ^P); r*l,2,... (3) 

with the boundary conditions: 

T^tP) *0 ? i=l,2,...,p; r*l,2,... (4) 

where A r is the rth eigenvalue and <j> I?) is the eigenfunction 
(sometimes also known as Node Shape; associated with X r 
Suppose the operator L is self-adjoint and positive definite, and 
all eigenvalues are positive and are ordered so that X^<^< ... . 
Since L is self -adjoint, the eigenfunctions are orthogonal and 
therefore can be normalized such that: 

/„>♦ (5) 

and + A,6fs (6) 

where 6^ is the Kronecker Delta. 

Using the esgacsicn theorem [3], the solution of u(P,t) can be 
obtained as: 

oo 

u(P,t) •JV'ttOiUt) O) 

r =1 

where 1^(0 is the nodal coordinate. Substituting (7) into (1), 
nultiplying both sides of the resulting expression by 
integrating over D and erploying (5) and (6), we obtain 

i^(t) Ujtt) ■ f r (t); r-1,2,... (8) 
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In (8) , the node (or natural) frequency w f in defined an 

u> r - b r ) V2 (« 

and the nodal control force f f (t) is ocaputed by: 

f r (t) - j^ff) f(p,t) ® (10) 

In practice, the infinite series in (7) in tnncated as 

u(P,t) « r (P) u^t) (11) 

where n is chosen to be sufficiently large so that u(P,t) can be 
represented with good fidelity. In this case we are dealing only 
with the first n nodes. 

Eg. (8) can be transformed into state equation form as 
follows: 


given below: 


ORIGINAL page » 
D£ POOR QUALITY 


" v ;] • 


r»l,2. 


( 25 ) 


Sines W f depends on ^ alone as seen in (20) , J can be 
arfniniied by ninlniring each J f independently. Frta optinal 
control theory (13) , the optinal solution for is given by 


C^(t) - R/ 1 ;r»l,2,....n (26) 

where ^(t) is the solution of the Riocati equation: 

i (t) « -* A - A T K ♦ K A' - Q ;r»l,2 n (27) 

r r r r r r r r 


x(t) = A x(t) ♦ wet) 

(12) 

From [4] the solution for G r (t) was obtained by 

x(t) - txj(t) x*(t)... xT(t))T 

(13) 

'■ ' 0 :j 

(28) 

W(t) = twj(t) wj(t)... tj(uf 

(14) 

ty b ♦ b) + r 1 


A * Block diag(A^, A^, ...,A r ) 

(15) 

r*- r/«* ) 

• * r p 

, (29) 

x r (t) - (u r (t) u r (t)/ Wr ) T 

(16) 

b - (i^. 2 + r f x ) l/2 ; r-1.2 n. 


W (t) « (0 1 

r~ n 

(17) 

We note that IMSC requires that the nntoer of actuators be 


and 


0 u> r 
~d) r o 


(18) 


for r=l,2,...,n. 

The system described by (12) consists of n subsystems given 

* 


x^t) * Aj. Xj.t + W r (t) ; r * 1,2,..., n. 


(19) 


The essence of IMSC is to choose V r (t) such that it depends 
on x^t) alone. Thus 


Wj(t) * Gj. r=l,2, ...,n 

where G r* prU «rl2l 
[ 5 r21 fl r22J 


r=l,2,...,n 


(20) 

( 21 ) 


equal to that of modelled modes. 

3. i ctumcr jmjammnm or Bnrsmrr kxal sna coma, 

Equations (19) and (20) represent the ooocept of IMSC. In 
order to implement (20), the modal state vectors x for 
r*l,2,...,n must be available for measurement. In [5] a r modal 
filter was developed to provide an estimate of the modal state 
vectors. Since this paper focuses on the problem of the 
implementation of actuators, we assume that the modal state 
vectors x r (t) are available for the state feedback in order to 
avoid the ocnplexity of getting involved in the state estimation 
problem that can well be the subject of a subsequent paper. 

Since it is impossible to control force at every point in 
the domain D, the distributed control force is realized by m 
(*<n) discrete point force actuators applied at m points P 2 , 
... ,P B in the domain D as given below: 


are (2x2) gain matrices. 

Substituting (16) and (17) into (20), we find that G r must 
asswe the following form: 


; r=l,2,...,n (22) 


For optimal control, g ^ and 0^22 ** determined such that 

the following quadratic cost function is minimized (linear 
regulator problem): 


ftp.t) ■ r t (t) 


(30) 


where 6(P-P«) is a spatial Dirac Delta fraction and F. (t) is the 
force applied by the ith actuator cn the point P^ . 1 

Now substituting (30) into (10) yields 


0 

0 

,®r21 

®r22. 


f (t) « /$ (P)L 6 (W.)F (t) dD 
r n r , , i 


(31) 


D r 1=1 

From the property of the Dirac Delta function, (31) can be 
reduced to 

m 


where 


faJflx + W T RV)dt 
r r r r r r 


(23) 

(24) 


Q r and IJ. are positive semidefinite and positive definite 
weighting matrix, respectively, associated with the rth mode. 

The form of G r given by (22) requires that \ assume the form 


\ t-(t) - Z p,(t). 

1=1 1 1 

If we define a force vector Fit) such that 

Fit) « t r x (t) r 2 (t) ... yt )] T , 

then using (32) the relatice between Fit) and V(t) can be 
expressed by 


(32) 


(33) 


W(t) - B F(t) 


(34) 
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ftm 8- |» 1 T *2 T ... B^ jT 

(35) 

|«(t)|p-||w-fif||^ (S-«f)Tsortf) 

s s 

(50) 

®Ui-l) * ®1»' r * 1 ' 2 " 

(36) 

** 91 *!(() * (sVrVsBw. 

(51) 

- » i (P l )/a i iiCty/eoi ... 

(37) 


(52) 


for r-l,2,...,n. is called the generalized inverse of F. 

4. FRCBLEH STATH®fT Proof. A proof of the above Lemma can be found in [12] • 

The implementation of the DISC scheme is illustrated in Figure 1 The main result of this paper is given in the folioring 

where the optimal state feedback law is defined by theorem: 


H(t) * G x(t) (38) 

•here 8<t) «= £Hj(t) Hpt) ... H^(t> ] T (39) 

and G * Block diagCGj, G 2 , G ). (40) 

The optimal solution for the DISC scheme was obtained in 
Section 3 as: 

V(t) * G x(t) . (41) 

The optimal solution is achieved if (41) is satisfied. In order 
to make V(t) equal to Gx(t) , the matrix 0 in Fig. 1 is designed 
such that W(t) « H(t) . From Fig. 1 we also have 

F(t) * D H(t) . (42) 

Kow substituting (42) into (34) , we obtain 

V(t) « BDKit). (43) 

In (43) to make W(t) * H(t), it is obvious that D is designed 
such that 

(44) 


From the structure of B as given by (35)-(37), each (2i-l)th 
row (odd row) of BD, for i*l,2,...,n is a row of zeros. We 
realise that (44) can never be satisfied. However r noting that 
each odd row of W(U is also a row of zeros, if we define 


B(t) *= (B* bJ... 

ir 

(45) 

6(t) * [d 2 d 4 ... 

■y 

(46) 


where Bj is the ith row of B and ^ the ith ooluon of D, then 
choosing a matrix D such that 

BD * L (47) 

will ensure that V(t) * Hit). It is noted that if (47) holds, 
then BD is a modified identity matrix of order (2nx2n) whose main 
diagonal elements are 0 at the (2i-l,2i-l) position and 1 at the 
(2i, 2i) position for i*l,2,...,n. 

One obvious solution for (47) is to choose D such that 
B * B . _ _ (48) 

However for the inverse of B to exist, B must be a nonsingular 
square matrix, requiring that the number of actuators be equal to 
the number of modelled modes («Fn) . Since the number of modelled 
modes is usually very huge for a faithful modelling of LSS, the 
required number of actuators is practically unrealizable. 

The problem considered in this paper can be formulated as to 
design the matrix D for a nonsquare matrix B (m<n) such that (47) 
is satisfied as well as possible. In other words, if this can be 
done, then the LSS can be controlled by a reduced number of 
actuators such that the closed -loop system is as optimal as 
possible. 

5. MUX RESULT 

bemmal: Consider the following equation: 

JMt) = BF(t) (49) 

where V(t) and F(t) are matrices consisting of even rows of W(t) 
and Fit ) , respectively. If the (nxa) matrix B has rank m, then 
the solution for (49) that minimizes the weighted norm of error 


Theorem 1 : Consider a large flexible space structure whose 
description and solution are given by (1) and (11) , respectively. 
If the operator L is self-adjoint, then there exists a control 
scheme with m (Kn) actuators that is suboptimal with respect to 
(24) in the sense that the closed-loop eigenvalues are assigned 
as closed as possible to those optimal eigenvalues specified by 
DOC. 

Proof: A control scheme with a recuced munber of actuators mild 
be optimal if B could be selected to be a ritfjt inverse of B in 
(47). However based on the form of B given ia (45) and (37) we 
can assume that 1 has rank m since the discrete actuators apply 
point forces at m distinct points Pj, P^-.-.P^- Fran [12] it is 
well-known that an (nxm) matrix (m<n) naving rank m does not 
possess any right inverse. Consequently B does not have any 
right inverse. According to Leona 1, because F(t) as given in 
(51) minimizes (50) , selecting a matrix 5 * B will minimize the 
difference jKHL, A making (47) be satisfied as sell as possible. 
Selecting 6 * 6* also make the closed-loop eigenvalues as 
identical as possible to those specified by DISC. Thus there 
exists a control scheme with a reduced ninber of actuators that 
is suboptimal with respect to (24). Q.E.D. 

6. EXAMPLE 

To illustrate the proposed control scheme we consider the control 
of a sinply supported beam whose dynamic is given by the Buler- 
Beraculli partial differential equation: 


U(dhx')ulx f t) + m(9 fet 2 )u(x,t) *f(x,t) (53) 

where for simplicity we set the mass m, the moment of inertia I, 
thae modulus of elasticity E and the length of the bean to unity. 
The boundary conditions fo this simply supported bean are: 
u(0,t) * u(l,t) * 0 (54) 

9 2 /3 ^u(0,t) * 3 ?3 x\i(l,t) * 0 (55) 

The solutions for the eigenvalue problem are given by: 

A j. « (Jdt) 2 (56) 

and w * TZ sin (kipc) . (57) 


Computer simulation was done to test the performance of the 
proposed control scheme. Results shew that the dynamics of the 
closed-loop system is affected considerably by the number and 
placement of actuators. Suppose we consider 20 modelled modes 
and divide the whole length of the beam into 20 sections 
specified by x(k)*k/21 for k-1,2,.,.,21. As expected, simulation 
results show that in every case the number of stable closed-loop 
eigenvalues increases with an increase in actuator number. The 
maximum numbers of stable closed loop eigenvalues when using 1, 
2, 3, 4 actuators are 5, 18, 19 and 20, respectively. The 
remaining eigenvalues are pure imaginary complex conjugate pairs, 
thus causing no instability but unwanted oscillations. It is 
also found that the actuator locations maximizing the number of 
closed loop eigenvalues are centered around both ends of the 
beam, namely at x*l/21 an & x*20/21. Fig. 2 presents the case for 
one actuator. As we see, the number of stable eigenvalues falls 
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off as the actuator acres toward the center of the ban. Similar 
results vara obtained for the cases of 2, 3 and 4 actuators. 
Table 1 prorides the dosed loop eigenvalues for DEC with 20 
actuators and the proposed control scheme with 4 actuators. As 
the table shows even with a reduced number of actuators the 
pr o posed control scheme assigns eigenvalues that are very close 
to those specified by the optimal DEC. Figures 3 and 4 shew the 
ben movement for the case of DISC with 20 actuators and the 
proposed control scheme with 4 actuators, respectively, when the 
ben is excited by an impulse. As the figures show, the proposed 
control scheme with a reduced ninber of actuators performs as 
well as the DEC with 20 actuators, bringing the bean movement 
down to zero after 5 second s . The time responses of mode 1 and 3 
are given in Figures 5 nd 6, respectively. We note that the 
time response of mode 1 of the control scheme is almost the same 
as that of DEC. There are seme insignificant differences in the 
time response of mode 2. The maximum difference is about 3x10 . 


7. (XMJJJSICH: 

In this paper, we first summarized the theory of DEC in the 
context of optimal control. Then the principal limitation of 
IMSC was pointed out in terms of the required number of 
actuators, being equal to that of modelled modes. After that, we 
proposed a control scheme which employs a reduced number of 
actuators to control a large space structure that in self- 
adjoint. Ccnputer simulation shewed that the number and placement 
of actuators play an important role in the stability of the 
closed-loop system. The research of this paper can be extended 
into studying the optimization problem of actuator placement and 
designing a weighting matrix S that maximizes the number of 
stable closed-loop eigenvalues. The problem of arbitrary 
assignment of eigenvalues [8] can also be addressed for the 
control of LSS. 
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Fig . 1 : Implementation of IMSC 



Fig. 2 : Relationship between the number 

of stable closed-loop modes and 
and the number of actuators 
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